Fermi- surface pockets in magnetic underdoped cuprates from first principles 
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Using an innovative first-principles band theory enabling the exploration of Mott-insulating mag- 
netic cuprates, we study the Fermi surface of underdoped Yi_ :E Ca a: Ba2Cu306 in a selection of mag- 
netically ordered and polaronic states. All phases exhibit qualitatively similar, hole-like nodal-point 
small pockets. Their properties (area, masses, mass sign) only partially match those extracted from 
recent quantum-oscillation experiments. Ab initio calculations, therefore, do not straightforwardly 
support a magnetic origin of quantum oscillations. 
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The existence of a Fermi surface (FS) in underdoped 
high-T c superconductors is attested pQ by angle-resolved 
photoemission (ARPES) [3J [5] and by recent Shubnikov- 
de Haas (SdH) and de Haas- van Alphen (dHvA) quantum 
oscillation observed in mildly underdoped (~0.1 holes 
per Cu0 2 unit) YBa 2 Cu 3 6 . 5 and YBa 2 Cu 4 8 [1 [5]. 
However, the FS detailed shape and origin continue to 
escape an unambiguous identification. Not only does 
SdH/dHvA indicate the existence of small Fermi pock- 
ets, whereas ARPES observes large zone-corner-centered 
cylinders at optimal doping turning into disconnected 
arcs near nodal points upon underdoping: in addition, 
the SdH/dHvA Hall resistance sign suggests electron-like 
pockets, whereas all ARPES measurements consistently 
report hole-like arcs. 

Reconciling arcs, pockets, and their electron vs hole 
nature is highly non trivial |6J. Ab initio calculations 
[5j [7] on the non-magnetic metal phase (the zero-order 
approximation to the normal state) of YBa2Cu30@.5 
and YBa2Cu40g indicate no solid evidence of FS pock- 
ets. Pockets may originate from a symmetry-breaking 
FS reconstruction, as suggested for "1/8" compounds 
[8] . Several models of symmetry breaking have appeared 
recently, among which <i-density-wave order [1], field- 
induced long-range magnetic order [9], short-range mag- 
netic order [TO] HH |T5], and magnetic polarons [131 
but none of them has yet gained a general consensus. 
Magnetic correlations are popular candidates as anti- 
ferromagnetic (AF) correlations coexist with, or survive 
into, the superconducting (SC) phase over a wide doping 
range, according to various experimental probes [T51ll6j . 

The very existence of a Fermi surface -hence of 
quasiparticles- and the hints about the role of magnetism 
suggest a considerable scope for ab initio band- like quasi- 
particle calculations, accounting accurately for material- 
specific information, and for coupling to the lattice. A 
first-principles description of underdoped cuprates, how- 
ever, is very challenging. Standard density functional 
theory [T7] fails to describe the magnetic insulating Mott 
phase and predicts a non-magnetic metallic state at all 
dopings, in strident conflict with the observed magnetic 



correlations. Here we overcome this limitation using the 
pseudo-self-interaction correction (pSIC) technique [TS] 
(already applied satisfactorily to many correlated mate- 
rials, including cuprates [HI [2Qj |2T] ) to calculate the FS 
of the AF and magnetic polaron states in an underdoped 
cuprate. All the states we investigate exhibit nodal pock- 
ets. In some cases, the pocket area is not far from ex- 
periment; the value and sign of the cyclotron masses, 
however, prevent a conclusive identification, leading us 
to conclude that the present configurations are not vi- 
able candidates for the observed oscillations. 

Here we study Y 1 _ 2 .Ca :I .Ba2Cu 3 6 at hole doping 
h=x/2=0.125. For this system hole injection only in- 
volves CuC>2 planes, without the complicacies [3T] due to 
oxygen doping. Calculations are carried out with a plane 
waves basis and ultrasoft [23] pseudopotentials (energy 
cutoff 30 Ryd, 12x12x12 special k-point grids for den- 
sity of states calculations, llxllxll uniform grid for 
Fermi surface calculations). Ca doping is described by 
explicit Y substitution in 2x2x1 supercells. 

The pSIC approach successfully describes [21] the un- 
doped precursor YBa 2 Cu30e as an AF Mott insulator. 
The main features of the AF Mott phase can be rec- 
ognized in the closely similar orbital-resolved density of 
states (OR-DOS) at h=x/2=0.125 (top panels of FigfTJ: 
valence and conduction bands are a mix of spin-polarized 
Cu & x 2_ y 2 and unpolarized O p^p^ states, with the main 
optical transition (involving Cu02-planc orbitals) start- 
ing at 1.2 eV, and higher transitions around 3.5 eV 
into apical O p z and Cu d z 2 states. The Cu magnetic 
moments of 0.5 \Ib and in-plane AF coupling constant 
J~0.15 eV in YBa2Cu306 agree well with experiment. 

Figs. [T] and [2] report our calculated OR-DOS, band en- 
ergies and FS for both AF and metallic non-magnetic (i.e. 
Pauli-paramagnetic (PM)) phases at h=0.125. The lat- 
ter shows doubly degenerate bands, and no Mott gap. Its 
valence band is a 3.5 eV-wide spin-unpolarized d x 2_ y 2- 
(Px,Py) hybrid, with Fermi level Ep at ~2 eV below 
the valence band top (VBT). The corresponding FS is 
a large cylinder centered at the zone corner, as expected 
for the metallic Fermi-liquid-like state. Only in-plane 
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FIG. 1: (Color on-line) OR-DOS for AF (top) and PM (bot- 
tom) phase at h=0.125 doping. Labels "ap" and "pi" indicate 
apical and in-plane (CUO2) atoms, respectively; arrows show 
spin directions. For each atom the light shadowed (green) 
curves shows the most important orbital contribution (p^ for 
O p i, p z for O ap , d x 2_ y 2 for Cu p i, Ba, and Y, d z 2 for Cu p z), the 
dark shadowed (black) curves the contribution of the other p 
(for O) and d (for Cu, Y, and Ba) orbitals. 



atoms contribute to the FS, whereas apical-atom states 
start appearing ^0.2 eV below Ep. In the AF state holes 
appear in the valence states of spin-polarized Cu, re- 
ducing their moments and weakening -though not yet 
destroying- AF order. The Mott gap splits the valence 
manifold into ~2 eV-wide empty upper and nearly filled 
lower Hubbard bands. The latter is sliced by Ep at 0.1 
eV below its top, along the BZ diagonal near the nodal 
point, with off-plane states again far below Ep. AF order 
causes a valence bandwidth collapse compared to the PM 
case, consequential to hindering nearest-neighbor Cu-Cu 
electron hopping. We remark that this result is a token 
of self-interaction removal in our calculations, and would 
not be found in a local-density- functional calculation 20J . 

A reason to compare the AF-phase FS with experiment 
is the evidence 24J that magnetic order develops in high 
fields in the superconducting state (this was modeled in 
Ref. [§])• The FS of the AF phase consists of small-area 
pockets around nodal points: in fact, four circlets appear 
around each node, one for each of the two up- and two 
down-polarized states per 2x2 CuC>2 plane. The average 
area per pocket converts to a frequency of 600 T, similar 
to the SdPT result in YBa 2 Cu 3 6 . 5 at h~0.1. The hole- 
like effective mass per pocket, however, is -0.5 m e , well 
below the experimental value. The carriers/holes ratio in 
the AF state at h=0.125 is n/h^l.2, suggesting a weak 
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FIG. 2: (Color on-line) Top: Band energies at x—1/4 
(h=0.125) for (lxl) PM phase (al) and (2x2) AF phase 
(bl). K-points coordinates (units of l/a=l/b, 1/c with 
a, b, c unit-cell parameters) are X=[7r/2,0,0], X=[0,7r/2,0], 
L=[7r/2,7r/2,0], M=[7r/2,7r/2,7r/2]. Bottom: calculated FS for 
PM (a2) and AF (b2) phase. 



violation of Luttinger's sum rule (n/h=l), or more likely 
that the FS differs slightly from the Luttinger surface 
to which the sum rule applies |25j . Experimentally jl] a 
similar value was deduced, n/h~1.5. 

Next, we take a first step towards modeling polaron 
states and the attendant FS, motivated by recent studies 
[P5] on strongly lattice-coupled polarons (a key concept 
in the original search for superconductivity in cuprates). 
In Fig|4] (again at h=0.125 doping) we explore singlet 
polarons of the Zhang-Rice type (ZRS) [26] embedded 
in the CuO^-plane AF background. The nature of these 
states is apparent from hole spin densities: in Fig.3(a2) 
we have one ZRS per cell (dashed circle), in Fig.3(b2) 
two ZRS per cell aligned into [110]-oriented stripes. AF 
symmetry imposes oxygen demagnetization and prevents 
ZRS formation, so we freeze-in an O breathing distortion 
towards Cu in each ZRS-to-be CUO4 unit (notice that 
the ZRS formation does not occur without distortion). 

We observe that the ZRS " condenses" on the distorted 
CUO4 plaquctte: holes localize on oxygens, and the in- 
duced magnetization is indeed anti-aligned to that of Cu, 
producing a vanishing total magnetization on the dis- 
torted CUO4 unit, as in the Zhang-Rice model. The ZRS 
signature in Fig.3(al) is a single CuC>2 band being de- 
pleted (i.e. hole-filled) and lifted by about 0.7 eV over 
the valence bands related to the AF-ordered units. In the 
stripe configuration of Fig.3(bl), two ZRS bands are now 
raised by roughly the same amount above the valence AF 
background. All other bands related to non-ZRS Cu02 
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FIG. 3: (Color on-line) Top: Band energies and hole spin 
density isosurface (h/V=0.006 bohr~ 3 ) at h=0.125 doping for 
two different magnetic phases with ZRS included, both within 
(2x2)symmetry. (al) and (a2) refer to the one-ZRS per 2x2 
plane configuration, (bl) and (b2) for two-ZRS per 2x2 plane 
configuration, (c) Orbital-resolved DOS for the one-ZRS con- 
figuration (only planar Cu d x 2_ y 2-0 (p x ,p y ) state are shown 
for simplicity). Red (dotted) ellipses enclose the ZRS DOS 
contributions from injected and native antialigned holes (blue 
arrows) to give S=0 locally. 



units are well below Ep and unaffected by hole injection. 
Thus, in both cases, only the ZRS bands contribute to 
the FS. The OR-DOS in Fig.3(c), left panel, shows that 
each ZRS band involves all the four in-plane oxygens first- 
neighbors to the central Cu of the ZRS unit. This up- 
polarized hole couples to the native Cu down-polarized 
hole (the two are enclosed in dashed ellipses) to form the 
singlet. The DOS from non-ZRS Cu0 2 units in (Fig.3(c), 
right panel, is close to that of the AF phase in Figjl] with 
hardly any magnetic moment reduction on Cu. 

Figj4] shows the FS of the three examined magnetic 
states (AF, one-ZRS, two-ZRS) at different values of h. 
The calculated FS at h=0.125 is in the central-column 
panels (the others are discussed below) . Closed FS struc- 
tures near nodal points appear to be a general feature of 
spin-polarized states. The one-ZRS FS shows a large 
nodal-point centered pocket due to the ZRS band (plus 
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FIG. 4: (Color on-line) Evolution of the FS vs. rigid Fermi 
level shift AEf for AF (top panels), AF plus one-ZRS (mid- 
dle panels), AF plus two-ZRS (bottom panels) configurations. 
From left to right, AE F =-0.2 eV (h~0.2), (h=0.125), and 
+0.1 eV (for AF), +0.5 eV for the two ZRS configurations 
(h~0.02). 

a minor AF-background pocket). In the two-ZRS con- 
figuration we find intersecting pockets from the two ZRS 
bands. The pocket shape is circular for the single ZRS, 
while for the stripe it is the intersection of two 90°- 
rotated ellipses. From the calculated FS areas of the ZRS 
bands, we obtain n/h~1.2 and ~1 respectively, again 
mildly deviating from the Luttinger sum rule. (In the PM 
and AF states, holes are evenly distributed on all CuC>2 
units. In the one- and two-ZRS configurations, they con- 
centrate on one and two Cu02 units respectively, so ef- 
fectively doping h=0.5 and h=0.25.) The frequency for 
the one (two) ZRS case is ^2500 (1600) T, with masses 
of -0.4 and -1 respectively. These values are not to close 
to experiment, and the mass sign is again hole-like. 

A definite conclusion we can draw is that different 
magnetic configurations will produce FS differing in the 
details, but sharing the character of nodal hole pock- 
ets. This brings us again to the fact that our calcu- 
lated FS pockets are hole-like, which is consistent with 
ARPES, but not with SdH/dHvA. To verify if electron- 
like branches may appear in the vicinity of Ep as calcu- 
lated for h=0.125, in Fig|4]we also plot the FS obtained 
rigidly shifting Ep, which roughly simulates doping fluc- 
tuations [7]. The FS remains hole-like and the pocket 
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area falls sharply for an upwards shift ("decreasing dop- 
ing", right panel column in Fig. 4). Electron pockets 
at the anti-nodes do instead appear (while hole pock- 
ets open up) for Ep shifting down (left panel column in 
Fig. 4). The shift, however, would correspond to opti- 
mal doping, i.e. to a unrealistic 50 % doping fluctuation. 
There is no obvious way to reconcile the Hall resistance 
sign with the present results. If one accepts the notion 
that the low T-high field Hall resistance (i.e. mass sign) 
change is a genuine electron current signature in the nor- 
mal state, all the present configurations should be ruled 
out. 

We finally touch upon attempts to reconcile SdH pock- 
ets -whose shape, number, and location is undetermined- 
with ARPES arcs near nodes (see e.g. [IT] Q2J [27j f° r a 
summary). Recently, a photoemission pattern switch- 
ing from large-doping Fermi cylinders to low-doping 
arcs was observed [5] upon tuning the self-doping of 
YBa2Cu306+rc via surface treatment. Cylinders and arcs 
appear to coexist in different sample portions depend- 
ing on the local self-doping, and arcs were hypothesized 
[31 H2] to be residual pockets (the estimated " virtual" ar- 
eas would be, as our own are, in fair agreement with SdH 
data). From total energies, we find the PM phase slightly 
lower than, but still quite close to AF. A coexistence 
may be expected, as suggested by /iSR experiments [15] . 
One may then speculate that the observed FS results 
from a configuration-averaged superposition of pockets 
and cylinder portions whose shape and spectral intensity 
depends on extension and shape of coexisting AF and 
PM regions. (Indeed, short-range AF order [TTJ [12] has 
been shown to cause a similar " unfocusing" of the pocket 
shape into an apparent arc shape.) 

In summary, we studied in detail the electronic struc- 
ture of underdoped Yi_ x Ca K Ba2Cu306 by an innovative 
first-principles band theory givin access to doped Mott 
magnetic states. AF phases (with or without embedded 
ZRS) produce small closed FS branches, while Fermi- 
liquid cylinders are related to spinless states. At low 
doping the spin-polarized holes are pervasive, and the FS 
is composed of pockets. The pocket/arc dichotomy may 
be due as a superposition of pockets and cylinder por- 
tions related to competing AF and PM regions. Some 
features of the pockets are consistent with oscillation ex- 
periments, but at least as many are not, so that no un- 
ambiguous identification can be made. We conclude that 
the doped antiferromagnet and the high-density ZRS po- 
laron configurations explored here are not those giving 
rise to the pockets inferred from experiments. It remains 
possible that more complex polaronic structures |13j be 
suitable candidates. 
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